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Dopant-based quantum computing implementations often require the dopants to be situated close 
to an interface to facilitate qubit manipulation with local gates. Interfaces not only modify the 
energies of the bound states but also affect their symmetry. Making use of the successful effective 
mass theory we study the energy spectra of acceptors in Si or Ge taking into account the quantum 
confinement, the dielectric mismatch and the central cell effects. The presence of an interface puts 
constraints to the allowed symmetries and lead to the splitting of the ground state in two Kramers 
doublets [J. Mol et al, App. Phys. Lett. 106, 203110 (2015)]. Inversion symmetry breaking also 
implies parity mixing which affects the allowed optical transitions. Consequences for acceptor qubits 
are discussed. 


I. INTRODUCTION 

Dopants, carrier providers for traditional transistors, 
are acquiring a more important role as semiconduc¬ 
tor nanostructures shrink their size. This change of 
paradigm was envisioned a few years ago with the first 
proposal for the implementation of a scalable dopant 
based quantum computer [T] . Silicon is an excellent plat¬ 
form for this technology due to its resilience against de¬ 
coherence and the already existing high level control of 
Si nanoelectronics [5]. Proposals relied originally on the 
quantum control of single electrons bound to donors, be¬ 
ing the electron spin [3], electron charge [4] or the nu¬ 
clear spin [T] the degrees of freedom used as qubits. Re¬ 
cent years have witnessed the practical demonstration of 
some of the ingredients involved in these proposals in¬ 
cluding single-shot spin readout of bound electrons [5], 
exchange coupling between electrons bound to neigh¬ 
bouring donors HE] and the practical implementation 
of the electrical control of spin qubits by an A-gate [S] as 
proposed in Ref. [1]. 

Proposals of acceptor-based qubits may make use of 
the long range strong dipolar inter-qubit coupling m , or 
exploit the spin-orbit interaction (which is stronger than 
for electrons) to couple spin to phonons [10], or to os¬ 
cillating electric fields m- The relative importance of 
the different sources of decoherence is different in elec¬ 
trons and holes: spin-orbit interactions [m dg would 
be more important for holes while hyperfine interaction 
(which can cause spin decoherence due to coupling to nu¬ 
clear spins [H] ) is smaller for holes than for electrons m- 
The effective suppression of the latter by Si isotopic pu¬ 
rification, which gets rid of nuclear spins in Si, leads to 
very long electron coherence times m- 

In practical dopant-based quantum computer propos¬ 
als, dopants are often introduced in nanostructures and 
close to surfaces or interfaces with materials different 
from the host. In this case, the energies of bound car¬ 
riers can be shifted by quantum confinement and dielec¬ 
tric mismatch [IlHISI potentially modifying the working 
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FIG. 1. Left: Schematic view of the problem. The acceptor 
A is at a distance d from the (001) interface between the 
semiconductor (Si or Ge) and an insulating barrier. Image 
charges appear due to the dielectric mismatch between the 
semiconductor and the barrier. Right: Sketch of the valence 
subbands. HH indicates the heavy hole subband while LH 
indicates the light hole subband. The energies of the bound 
states are positive and defined with respect to the top of the 
valence band. 


parameters of the devices. Quantum confinement may 
alter the shape of the wave-functions through the bound¬ 
ary conditions, consequently affecting the binding energy. 
For instance, in a very thin (compared to the bound state 
wave-function size as) nanowire, the extra confinement 
enhances the binding energy deactivating the dopants as 
carrier providers uniiiii- However, when the dopant is 
close (compared to as) to one interface/surface but not 
confined in other directions, the wave-function can be 
deformed in such a way that the density probability of 
the bound state decreases on the dopant, leading to a re¬ 
duction of the binding energy. The dielectric mismatch 
gives rise to image charges which, depending on the rel¬ 
ative size of the dielectric functions of the nanostructure 
components, can lead to an enhanced or decreased bind¬ 
ing. In the case of a semiconductor surrounded by in¬ 
sulators, the image charges have the same sign as the 
charges originating them, enhancing the binding ener¬ 
gies. Not only the energies but also the symmetry of 
the bound states may be modified. In the case of accep¬ 
tors, the four-fold degeneracy of the ground state may 
be broken by strain electric field [23 [22] and/or 
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magnetic fields. Characterisation of dopants embedded 
in nanostructures can be performed via transport mea¬ 
surements miiTiiiHi while Scanning Tunnel Microscopy 
can give information of the wave-function of subsurface 
dopants [2M52] . 

Here we perform an analysis of the energy spectrum 
and the symmetry of the bound states for substitutional 
acceptors (group III elements) in Si or Ge close to an 
interface with an insulating barrier. Our study is based 
on effective mass theory (EMT). EMT, despite its ap¬ 
parent simplicity, has been proven to give very accurate 
descriptions not only of the binding energies but also of 
the donor wave functions [32] . EMT exploits the analogy 
with free atoms but includes information about the host 
crystal through the bands. Both quantum confinement 
and dielectric mismatch are included. Central cell cor¬ 
rections |33l [34] are also considered in order to reproduce 
the energy spectra of different acceptor species. A recent 
publication |3T] has demonstrated that the four-fold de¬ 
generacy of the ground state is broken for B acceptors 
in close proximity (~ nm) to a surface. We study the 
energy spectra of acceptors in Si and Ge, including the 
observed Kramers doublet splitting of the ground state. 
Our approach allows a complete analysis of the symme¬ 
try breaking induced by the interface. Besides a symme¬ 
try reduction analogous to the one produced by uniaxial 
strain |23j , the interface also breaks inversion symmetry, 
leading to parity mixing. The consequences for an accep¬ 
tor based qubit are analysed. 


II. MODEL 

The energy levels of a single acceptor are mainly deter¬ 
mined by the valence band structure of the group IV host 


crystal, see Fig. The valence bands have their max¬ 
imum at the r point and the main atomic contribution 
belongs to p states. This implies that the total angular 
momentum of the bands can be J = 3/2 or J = 1/2. 
J = ?>12 corresponds to the T^ bands. At fc = 0 these 
bands are four-fold degenerate but split away from the T 
point into two doubly degenerate bands: The heavy-hole 
bands (with |mj| =3/2) and the light-hole bands (with 
|mj| = 1/2). The T^ band corresponds to J = 1/2 and 
it splits from the T^ bands by the spin-orbit coupling 
^SO- 

The Kohn-Luttinger Hamiltonian [3^ describes a posi¬ 
tive hole bound to an acceptor centre taking into account 
the six valence subbands. The Hamiltonian is a 6 x 6 
matrix operator which accounts for all the possible com¬ 
binations of the quantum numbers J and Jz- In Si, the 
spin-orbit coupling is Aso = 44 meV |36] which is com¬ 
parable to the energy of the acceptor ground state (> 45 
meV), and hence the contribution of the split-off band 
can not be neglected. This is not the case for Ge, where 
the split-off energy is 290 meV [33] , much larger than the 
binding energy (~ 10 meV), allowing in principle the re¬ 
duction of the dimensionality of the matrix operator to a 
4x4 matrix. In the following we will keep the full 6 x 6 
Hamiltonian for both Si and Ge for completeness. 

In bulk, it is possible to reformulate the Kohn- 
Luttinger Hamiltonian separating the spherical symmet¬ 
ric terms, like the Goulomb impurity potential, from the 
terms with the cubic symmetry of the crystal [37l439j . 
This allows an easier analysis of the symmetries and se¬ 
lection rules of the acceptor states. However, this simpli¬ 
fication is not possible in the presence of an interface due 


to the reduction of the symmetry, as explained in Sec. HI 


Hence we keep the original 6 x 6 Kohn-Luttinger Hamil¬ 
tonian [33] 
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Dehning the effective Rydberg unit as Ry* 
e^mo/ 2 Ii^e^ 7 i and the effective Bohr radius as a* 
h'^Csji/e^mo, the differential operators in Eq. Q are 

P=-e+- 


Q = -ll{kl + kl-2kl) 

7i 


L = i2'/3—{kx — iky)kz 
7i 

M = -V3^ikl-kl)+i2V3^k,ky, (2) 
7i 7i 

with Too the free electron mass, £„ the semiconductor 
static dielectric constant, and 71 , 72 and 73 material de¬ 
pendent Luttinger parameters [35] related to the curva¬ 
ture of the valence subbands. For Si Ay* (Si) = 24.8 
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meV, a*(Si) = 2.55 nm, 7 i(Si) = 4.27, 72 ( 8 !) = 0.32, 
73 (Si) = 1.458, while for Ge Ry*{Ge) = 4.4 meV, 
a*(Ge) = 10.85 nm, 7 i(Ge) = 13.35, 72 (Ge) = 4.25, 
73 (Ge) = 5.69 [40] . 

In order to account for the dependence on the accep¬ 
tor species of the binding energies, the so-called central 
cell corrections have to be included [33] . We adopt here a 
central cell potential which takes into account the incom¬ 
plete screening of the Coulomb potential at very short 
distances from the dopant [34] 


Kc = 



(3) 


with Tec a semiempirical parameter [34] calculated such 
that for a given acceptor the measured bulk ground state 
energy is reproduced, see Table jl] A single Tcc charac¬ 
teristic of each dopant species is sufficient to get also 
the excited spectrum. The values are very similar for 
Si and Ge. However, due to the much smaller binding 
energies of the acceptor states in Ge (which corresponds 
to much more extended wave functions) the effect of the 
central cell correction on the binding energies is not as 
large in Ge as in Si. The central cell correction is not 
needed to reproduce the energy spectrum for the boron 
acceptor (namely, Tcc for boron is negligibly small). A 
larger binding energy corresponds to a larger Ccc- Typi¬ 
cal central-cell parameter values are very small, Tcc 1 
A, and hence we do not expect it to be affected by the 
presence of the interface. 




B 

A1 

Ga 

In 

T1 


Egs (meV) 

45.83 

69.03 

74.16 

157 

246 

Si 


rcc (nm) 

- 

0.078 

0.082 

0.12 

0.15 


Egs (meV) 

10.82 

11.15 

11.32 

11.99 

13.45 

Ge 


rcc (nm) 

- 

0.077 

0.089 

0.12 

0.16 


TABLE I. Central cell parameter Tcc |33| that reproduces the 
measured bulk ground state energy Eqs for the different ac¬ 
ceptor species [3S]. A single suffices to reproduce the full 
spectrum [34|. Boron binding energies are well reproduced 
without the central cell correction. 


The effect of the proximity to the interface is consid¬ 
ered by including in the Hamiltonian the image charges 
that arise due to the dielectric mismatch between the 
host crystal and the barrier [HUH]- A (001) interface is 
considered at a distance d from the acceptor, see Fig. [^ 
The total Hamiltonian is then 




acceptor 


= Hkl + Vec — 


2Q' 


+ 


Q' 


+ {z + 2d)2 2{z + d) 


(4) 

where the third and fourth terms are the acceptor and 
the hole images respectively with Q' = {ei, — es)/{eb + es). 
tb is the barrier static dielectric constant m- Note that 


for an insulating barrier Q' < 0, namely the acceptor 
image is attractive for holes, and hence an enhancement 
of the binding energy is expected. Typical barriers con¬ 
sidered have very large gaps compared to the binding 
energies involved. Therefore, we will assume the bound 
hole encounters a hard wall at the interface. We neglect 
the electric field potential in order to focus on the effect 
of the symmetry reduction by the interface. 

The Hamiltonian in Eq. Q has the symmetry of the 
bulk crystal. Close to an interface, the symmetry of the 
crystal is reduced. Effective changes in the spectrum, 
similar to the ones caused by uniaxial strain with the 
splitting of the light-hole and heavy-hole subbands [331 
1431144] . are expected. These are discussed in the next 
section. 


III. SYMMETRIES AND VARIATIONAL 
METHOD 


In bulk, the acceptor problem has cubic symmetry with 
inversion so the transformation elements form the Oh 
group. The three irreducible representations Tg and Fy 
with dimension two and Fg with dimension four are all 
allowed in this group. As a result, the acceptor states 
are doubly or four-fold degenerate. Due to the inversion 
symmetry, parity is conserved and the states can be sep¬ 
arated in well defined parity states (F)!", Fg ^ Fy , Fy , F^ 
and F^). However, central cell effects break the inver¬ 
sion symmetry, reducing the symmetry of the system to 
the tetrahedral double group T^. The T^i group also has 
irreducible representations Fg, Fy and Fg but in this case 
parity is not a good quantum number. A consequence 
of the inversion symmetry breaking is the appearance of 
weak transitions between states with nominally identical 
parity, which would be forbidden if inversion symmetry 
were actually preserved [38|. However, as the central cell 
corrections are only important at small distances from 
the acceptor (see Table |l]), their effect is very local and 
hence the parity can be in general considered a good 
quantum number. In the Hamiltonian Eq. Q, we ne¬ 
glect this symmetry breaking as we consider a central 
cell correction with spherical symmetry. 

The selection rules of the Hamiltonian Eq. Q can be 
obtained after examining the differential operators and 
the couplings between different sets of the pseudo-angular 
momentum | J, J^). Being the atomic orbital angular mo¬ 
mentum L, we can define a total angular momentum F as 
F = L -I- J [3H1139| . The selection rules for the quantum 
numbers F and which can be obtained ap¬ 

plying the Wigner-Eckart theorem to the terms of cubic 
symmetry, are 


|77acceptor 1^5 As) OC j . (5) 

We can use these selection rules to relate the quantum 
numbers F and to the irreducible representations of 
the Oh (or T group. Doubly degenerate eigenstates 
of the cubic symmetry with Fz = F.1/2 F An transform 
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under the group as Fg states, while two-fold degenerate 
states with = ±3/2 ± 4n transform like Fy symmetry 
states. The four-fold degenerate eigenstates of the cubic 
symmetric terms in the Hamiltonian correspond to the Fg 
representation and can have any half-integer F^, always 
according to the selection rules. 

The description of the states is modified in the presence 
of an interface. The inversion symmetry is broken by the 
interface, here assumed to be in the ( 001 ) direction, so 
the parity is clearly not a well defined quantum number. 
Moreover, in the semi-infinite space, the spherical har¬ 
monics do not form an orthogonal basis and therefore L 
can not be a well defined quantum number. An imme¬ 
diate consequence is that the total angular momentum 
F = L ± J is not well defined and states with differ¬ 
ent F are not orthogonal to each other. However, the 
z projection of the atomic angular momentum is asso¬ 
ciated to the (fi coordinate so Lz is not affected by the 
presence of the interface and, as Fz is the sum of Lz and 
Jz, the selection rule Eq. ([^ holds. In terms of symme¬ 
try, the tetrahedral double group is reduced to the 
tetragonal group D 2 d- The Fg symmetry is not allowed 
in the D 2 d group but the irreducible representations Fg 
and Fy remain. This implies that close to an interface 
the four-fold degeneracy of the Fg states is broken into 
two doubly degenerate states with symmetries Fg and Fy 
respectively. This effect of symmetry reduction by the 
interface is analogous to the effect of uniaxially strained 
silicon in the (001) direction [2^ [24l l45] . 

The barrier potentials usually considered are much 
larger than the typical binding energies and hence a hard- 
wall boundary condition 4'(z < — d) = 0 is assumed for 
the wave function. The interface boundary condition im¬ 
plies that the spherical symmetry usually assumed for 
the bound hole variational wave-function in bulk is not 
valid and it is more appropriate to work in cylindrical 
variables with the z-axis perpendicular to the interface. 
With the information of the symmetries and the selection 
rules we can define a hydrogen-like variational basis set 

V-(p, z, ¥>, a,) = (z ± J, Jz) , ( 6 ) 

where V = L—\Lz \ and n' = n — L—\ with n > L. A set 
of different Oi valu^is considered, p is the cylindrical 
variable p = \/r^ — . The (z ± d) prefactor assures that 

the basis set satishes the hard wall boundary condition. 
Given an Lz , the quantum number Jz is chosen such that 
Fz belongs to the convenient symmetry. This basis set is 
truncated at a certain Lmax with Umax = Lmax + ^- L^ax 
is chosen such that the condition for the energy — 

= 0.1 meV is fulfilled in bulk, which in silicon 
corresponds to Lmax = H- The number of different 
considered is not as determinant for the convergence as 
the value of Lmax- For example, for the ground state 
energy of an acceptor in bulk silicon, including oi = 1 
and 02 = 2 with Lmax = 11 gives - El^^^-i) = 

0.1 meV. However, the excited states require a larger set 
of different ai due to their different extensions. Adding 
og = 0.5 and 04 = 0.25 gives results with — 


< 0.2 meV for the first 8 states. For Ge, the 
set is truncated at Lmax = H for the J = 3j2 states 
while for the J = 1/2 states Lmax = 6 since the split-off 
band is less relevant in this case and the convergence is 
faster. 

Since the basis set defined in Eq. Q is not orthonormal 
it is necessary to consider the overlap matrix Sij 

^i,j — {'4^(^^i')\^acceptor\'^{^j'}) : (7) 

and the problem becomes a generalized eigenvalue prob¬ 
lem 

( 8 ) 

Most of the integrals used to obtain the matrix ele¬ 
ments of Hij and Sij can be solved formally, as detailed 
in Appendix [A] 

IV. RESULTS AND DISCUSSION 

Fig.@ illustrates the effect of the interface on the en¬ 
ergy spectrum as a function of the interface-acceptor dis¬ 
tance d in Si. The superimposed dots at large d corre¬ 
spond to the three lowest energies measured in bulk for 
the corresponding acceptors [31EO]. The bulk energies 
are well reproduced but the interface boundary condition 
and the image charges affect the calculated energies up to 
the distances shown in the hgure. Fig. (a) corresponds 
to a B acceptor and neglects the image charges (namely, 
Q' = 0), Fig. [^ (b) considers a Si 02 barrier {Q' = —1/2), 
and Fig. (cjis the result for Al acceptors with Q' = 0. 
Fig. [3] shows the corresponding results for B acceptors in 
Ge with Q' = 0. There are two main qualitative interface 
induced effects on the energy levels: (i) the binding ener¬ 
gies are smaller close to the interface and (ii) the ground 
state (which is four-fold degenerate in bulk) splits in two 
Kramers doublets |31| . 

The reduction of the binding energies close to the inter¬ 
face is due to the quantum confinement [IMS] produced 
by the boundary condition on the wave function, which 
has to be zero at the interface. The wave-function is 
hence deformed with its probability density shifting away 
from the interface, see Fig. This effect is more signifi¬ 
cant for the levels coming from the four-fold degenerate 
bulk ground-state than for the excited states leading to 
the compression of the full energy spectrum [46] . This 
compression also appears in bulk strained systems 
due to the splitting of the heavy-hole and the light-hole 
bands. 

The reduction of the ground state binding energy due 
to the quantum conhnement is partially compensated by 
the dielectric mismatch with the insulating barrier m, 
compare panels (a) and (b) in Fig. the holes are more 
strongly bound when the image charges are included (be¬ 
cause Q' < 0). The extra binding effect of the acceptor 
image can be still appreciated at the longest distances 
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FIG. 2. Energy dependence on distance to the interface of the ground state and some excited states of an acceptor in Si. The 
highest energy corresponds to the ground state which is four-fold degenerate in bulk and split in two Kramers doublets near 
the interface. Dashed lines indicate Fy states while solid lines indicate Fe states. The red dots are the experimental values for 
the bulk energies of the lowest three states [3S1I10]. As the excited states have a less localized wave function than the ground 
state, they become affected by the interface at larger distances, (a) For B acceptors, neglecting the image charges, (b) Same 
as (a) including the image charges corresponding to a Si02 barrier. Note the enhancement of the binding energies in (b) with 
respect to (a) due to the attractive character of the acceptor image, (c) Same as (a) for A1 acceptors. Holes are more strongly 
bound for A1 in bulk, see Tablebut the energies at small d are very similar to the B case in (a). 



FIG. 3. Spectrum of a B acceptor in Ge. The notation 
for the curves is the same as in Fig. The binding energy 
is much smaller in Ge than in Si, however, the ground state 
splittings are of the same order, see Fig. 

shown in Fig. (b) by comparing to the bulk values. A 
vacuum barrier, with Q' = —0.84, would further increase 
the ground state binding energy such that, for d = 2 nm, 
Egs = 40.3 meV, consistent with the reported experi¬ 
mental values in Ref. [19]. 

For doped Si, the energy difference between the two 
doublets and the excited spectrum can be lowered (for 
Q' = 0) to values < 8 meV, a significant reduction from 
the > 25 meV typical splitting in bulk. This energy 
splitting is enhanced when the dielectric mismatch is in¬ 
cluded, as in Fig. (b), but still smaller than its bulk 
value. The compression of the spectrum should be kept 
in mind when interpreting experimental measurements of 
bound states in field effect transistor geometries where 
Si 02 is a common barrier material. 

Qualitatively similar results are found for A1 acceptors 
in Fig. (c). In the bulk limit, the energies of the two 
first levels are enhanced with respect to B by the central 
cell corrections. However, the third level has a binding 
energy very similar to that in B. This difference is due to 
the first 2 energy levels being s-like (and hence more af¬ 
fected by central cell effects). Notoriously, although the 


A1 acceptors have much larger binding energies in bulk, 
close to the interface the values are very similar to B 
acceptors. This is a consequence of the hole probabil¬ 
ity density shifting away from the dopant, significantly 
reducing the effect of the central cell correction on the 
energy. Therefore, distinction among different acceptors 
in terms of measured binding energies may be blurred by 
the proximity to an interface. 

The splitting of the ground state is due to the symme¬ 
try reduction produced by the interface. As explained 
in Sec. |III[ the Fg symmetry is not allowed and hence 
the states with Fg symmetry in bulk acquire a Fg or r 7 
symmetry as d is reduced. This is the case in particu¬ 
lar of the four-fold degenerate bulk ground state which is 
split in two doublets. The two doublets have a predom¬ 
inant light-hole (Fg) or heavy-hole (r 7 ) character. Both 
corresponding wavefunctions have s-like envelopes, how¬ 
ever, the light-hole ground state is more affected by the 
interface because it has a higher contribution of high or¬ 
der spherical harmonics parallel to the surface as shown 
in Fig. [^ The slightly different shapes of the heavy-hole 
and light-hole wave functions leads to the energy splitting 
of the two doublets [5T] . 

The solid and dashed lines in Fig. correspond to Fg 
and r 7 symmetries respectively. Whenever those curves 
are degenerate towards the bulk (increasing d), the Fg 
symmetry is recovered. The level crossings (anticross¬ 
ings) in the excited spectrum occur between states with 
different (same) symmetry. For small values of d there are 
some near degeneracies between Fg and F 7 states which 
are accidental and not related to the (prohibited) Fg sym¬ 
metry. 

Fig. [5] shows the energy difference between the two 
lowest doublets for B in Si. The energies for this plot 
have been obtained considering an interface with vacuum 
{Q' = —0.84) but results are independent of the value of 
Q' considered. The energy difference is enhanced as the 
acceptor gets closer to the interface. A decreasing split- 
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FIG. 4. Spatial probability distribution of the two ground 
state doublets for B in Si considering a Si02 interface: (left) 
Fe state with a predominant light-hole (LH) character, and 
(right) Ft state with a predominant heavy-hole (HH) char¬ 
acter. The acceptor is located at (0,0,0). Top figures are 
the in-plane {x — y) images for d = 2 nm. The shape of the 
wave-functions is the same at all distances. Differences with 
d can be noticed in the x — z plane (images are equivalent in 
the y — z plane). The middle panels correspond to d = 7.5 
nm and the bottom ones to d = 1.5 nm. The red line in the 
bottom panels represents the interface position. The wave- 
functions are deformed by the interface proximity. The LH 
wave-function is more affected than the HH one leading to 
the energy splitting. 


ting is found for distances dmin ~ 0.5a*, smaller that the 
ones shown in the figures, probably signalling a qualita¬ 
tive change in symmetry similar to the one found in the 
limit of infinite strain, with extra degeneracies [2311211 ■ 
However, we do not expect the effective mass approxima¬ 
tion used here to be valid for values d < dmin- 

The dots with error bars in Fig. are experimental 



d(nm) 

FIG. 5. Splitting of the B ground state in Si as a function 
of the distance d. The red dots are the experimental values 
taken from Ref. m- Although an interface with vacuum has 
been considered in this plot {Q' = —0.84), the value of the 
splitting is basically independent of Q' down to the distance 
d considered here. 



FIG. 6. Comparison of the splitting of the ground state in 
Si (with two different acceptors) and Ge in effective units of 
distance. 


values as reported in Ref. m- The agreement with the 
calculated splitting is very good even in comparison to 
more sophisticated tight-binding calculations |3T]. The 
giant splitting ^ 6 meV found at very short distances 
d ^ 0.5 nm < dmin in Ref. |3I] is not included in the 
figure. We stress that the doublet splitting does not de- 



FIG. 7. This figure illustrates how the parity is lost in the 
proximity of an interface for B in Si and neglecting the image 
charges {Q' = 0). The color code for the lines is the same as 
in Fig. The solid lines indicate the Fe states while dashed 
lines indicate Fy. The black and orange lines are the first two 
even states (bulk ground state). Red and blue are the first 
two odd states. 
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pend on Q', namely it is independent of the nature of the 
insulator, but the energy spectrum compression does, see 
discussion of Fig. Therefore, in order to keep the ex¬ 
cited states away from the lowest doublet in case of a 
giant splitting, an insulator with a very low dielectric 
constant should be used. 

The same considerations can be made in the case of 
doped Ge, see Fig. for the energy spectrum. The main 
difference with Si is that holes are much less bound in 
Ge, leading to a larger effective Bohr radius and hence 
the effect of the interface is notorious for much deeper 
acceptors. 

Although the energy spectrum is quantitatively af¬ 
fected by the host crystal and the acceptor species, the 
interface induced ground state splitting is very similar 
in different systems, as illustrated in Fig. Here the 
doublet splitting is shown for B in Si, A1 in Si and B in 
Ge as a function of the distance in effective units. The 
three curves are very similar and in all cases splittings 
> 1 meV can be found. In comparison, the electric field 
required to achieve this splitting is much larger than the 
ionization field 5 MV/m justifying the omission of 
the electric field in our model. 

As explained in Sec. |IIH breaking the inversion symme¬ 
try implies that the parity of a state is not a conserved 
quantity anymore. Fig. quantifies (for the particular 
case of B in Si) how the parity of the first two even and 
two odd states lose their well defined parities. The wave- 
function tk is a combination of the basis states ^ in Eq. 
each with a well defined parity P = ±1. The parity in the 
figure is defined as ('I'|P|'I')/(vI/|v]/). Even at the longest 
distances shown, d = 25 nm, the different states are not 
completely even or odd. This implies that states with 
different parities have a finite overlap {S ^ 5% already 
at d = 25 nm). S increases strongly with decreasing d, 
reaching S > 0.2 at d ~ 2 nm. Consequently, optical 
transitions which are not allowed in bulk, may be per¬ 
mitted close to an interface. 


V. SUMMARY AND CONCLUSIONS 

We have used an effective mass approach to study the 
effect of an interface on the group III acceptor energy 
spectra and bound states symmetry in Si and Ge. This 
method, though computationally less demanding than 
tight-binding, has been proven to be very reliable not 
only for the calculation of binding energies but also for 
the wave-functions [32]. A semiconductor/insulator in¬ 
terface introduces a specific boundary condition and the 
corresponding dielectric mismatch implies new attractive 
potential terms in the Hamiltonian. The combined effect 
of the quantum confinement and the dielectric mismatch 
gives rise to energy shiftings which depend on the dis¬ 
tance between the acceptor and the interface and the 


value of the insulator dielectric constant. In general 
terms the full energy spectrum is compressed, namely, 
the distance between the ground state and the excited 
states is reduced by the confinement but this reduction 
is partially compensated by the dielectric mismatch m- 
Gentral cell effects, which account for the binding en¬ 
ergies dependance on the acceptor species, become less 
important when acceptors get closer to the interface with 
the insulator. Therefore, different acceptors close to an 
interface may be difficult to distinguish by the values of 
their binding energies. 

We have also followed the modifications on the sym¬ 
metry of the bound states which can be qualitatively 
understood via the analysis of the symmetry breaking 
induced by the interface. One of the consequences of 
this symmetry reduction is the splitting of the four-fold 
ground state in two Kramers doublets, as reported in 
Ref. [31]. This doublet splitting is independent of the 
dielectric mismatch for a particular host:acceptor com¬ 
bination. Different acceptors lead to comparable values 
of the interface induced doublet splitting, which is typi¬ 
cally > 1 meV. Our results are in very good agreement 
with the measurements in Ref. m- Inversion symmetry 
breaking implies that the parity is not a good quantum 
number leading to parity mixing and allowing for optical 
transitions which are not allwed in bulk. 

It will be useful to keep in mind these results in the 
implementation of an acceptor based quantum computer. 
The strong dependence of the binding energies, particu¬ 
larly for the ground state, on the distance between the 
acceptors and the interface introduces an uncertainty in 
the sample preparation due to the difficulty in control¬ 
ling the positioning of dopants to the required degree. 
The interface induced ground state splitting could any¬ 
way be exploited after a posteriori characterisation of 
the samples. Strain could provide a way to induce the 
ground state splitting if acceptors are buried farther from 
the interface to avoid acceptor to acceptor binding en¬ 
ergy variability. Strain would also be more advantageous 
than interface effects if optical transitions are involved, 
as the former does not produce the parity mixing we find 
for subsurface acceptors. In any case, the fact that the 
ground state splits due to interface proximity allows the 
possibility of a tunable two-level system in state-of-the- 
art devices. 
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Appendix A: Integrals 


The calculation of the matrix elements Hi j and Sij have required the evaluation of integrals involving products of 
wavefunctions (§ and the expected values of the different operators (§. The analytical solutions to these integrals 
are not tabulated. We summarize here most of the integrals used in this work. 

The general form of the integrals involved in the calculations of the matrix elements is 


/ oo poo _ 

dz J dr ^\/j 


k' 


(Al) 


The evaluation of this integral depends strongly on the parity of the exponent k'. When k' = 2k being k any 
positive integer and defining for convenience j = 2k + n + c + 2: 


I{c,n,2k) = 


k\T{^) 
2T{^ + k) 


(r(7)(i + (-ir)) + (-ir^r(7,ad) 


being r(a, z) the incomplete gamma function. 

When k' is an odd number k' = 2k + the parity of the exponent n becomes relevant. If n is even 
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being F the hypergeometric regularized function. And when n is odd 

"y + l — c 


i{c.n.F) = j2(-ir Y. dl^' + i) 

m—0 u—1 l—O ^ 
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(A3) 


(A4) 


The exponent n can be negative for certain operators. When this is the case, the integrals can be transformed into 
the previous integrals by using the method of differentiation on a under the integral sign. For example, when n = — 1 


/(c, —1, fc') = — y da I{c,0,k') + C, 


(A5) 


where (7 is a constant that can be obtained using I{c, n,k) = 0 when a —>■ 00 . The rest of negative n integrals can be 
obtained using this method recursively. 
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